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Abstract 

We consider the linear elliptic equation —div(aVu) = / on some bounded domain D, where 
a has the form a = exp( 6 ) with h a random function defined as h{y) = where 

y = iVj) S i.i.d. standard scalar Gaussian variables and is a given sequence of 

functions in L°°{D). We study the summability properties of Hermite-type expansions of the 
solution map y i—> u(j/) € C := Hq{D)^ that is, expansions of the form u{y) = u^H^^y), 

where Hu[y) = nj>i Hif^{yj) are the tensorized Hermite polynomials indexed by the set T of 
finitely supported sequences of nonnegative integers. Previous results [19] have demonstrated 
that, for any 0 < p < 1, the P summability of the sequence (jljV't ||L°°)i>i implies P summability 
of the sequence (||u,y||y),ygjr. Such results ensure convergence rates n”® with s = ^ ^ of 

polynomial approximations obtained by best n-term truncation of Hermite series, where the 
error is measured in the mean-square sense, that is, in L^(K^,y, 7 ), where 7 is the infinite¬ 
dimensional Gaussian measure. In this paper we considerably improve these results by providing 
sufficient conditions for the P summability of (||uj/||y),yg:F expressed in terms of the pointwise 
summability properties of the sequence i\tpj\)j>i- This leads to a refined analysis which takes 
into account the amount of overlap between the supports of the ipj. For instance, in the case of 
disjoint supports, our results imply that, for all 0 < p < 2 the summability of (||mi,|| 
follows from the weaker assumption that (||V'jlli°°)i>i is summable for q := > p. In the 

case of arbitrary supports, our results imply that the summability of (||u;/||y);/g:F follows from 
the (P summability of (j^jj'bjjj for some /3 > ^ which still represents an improvement over 

the condition in [19]. We also explore intermediate cases of functions with local yet overlapping 
supports, such as wavelet bases. One interesting observation following from our analysis is that 
for certain relevant examples, the use of the Karhunen-Loeve basis for the representation of b 
might be suboptimal compared to other representations, in terms of the resulting summability 
properties of (|[u,y|[y),yg:F. While we focus on the diffusion equation, our analysis applies to 
other type of linear PDEs with similar lognormal dependence in the coefficients. 
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1 Introduction 


1.1 Approximation of high-dimensional parametric PDEs 

Parametric partial differential equations equations have the general form 

P(u,y)=0, (1) 

where u i—)• T){u,y) is a partial differential operator that depends on a parameter vector y = 
ranging in a certain domain U G where J > 1 is hxed. Assuming well-posedness of 
the problem in some Banach space V, the solution map 

y ^ w(y), (2) 

is dehned from the parameter domain U to the solution space V. 

Equations of this type arise both in stochastic and deterministic modeling, depending on the 
nature of the parameters yj which may either be random or deterministic variables. In both settings, 
one main computational challenge is to approximate the entire solution map y e-)• u{y) up to a 
prescribed accuracy, with reasonable computational cost. This task has been intensively studied 
since the 1990s, see in particular [14, 15, 20, 22] for general treatments, ft becomes very challenging 
when the number of parameters J is large due to the curse of dimensionality. Ideally, one would 
like to design numerical methods that are immune to the growth of J, which in principle amounts 
to treat the case of countably many variables, that is, 

y = {yj)j>i G P c (3) 

This problem has been the object of much attention in recent years [3, 4, 7, 8, 9, 16, 13, 19]. 

Sparse polynomial methods are based on approximations to u of the form 

UK[y) ■='^Uyy'', (4) 

where A C T" is a hnite set of (multi-)indices v = {vj)j>i G T and y^ = OjM y'^/ ■ the case of 
an inhnite number of parameters, the index set F denotes the (countable) set of all sequences of 
nonnegative integers which are finitely supported (i.e. those sequences for which only finitely many 
terms are nonzero). Note that the polynomial coefficients are functions in V, and therefore the 
construction of ua requires in principle the computation of #(A) such functions. 

One particularly relevant example is the model elliptic PDE 

—div(aVrt) = /, (5) 

set on a bounded Lipschitz domain D <Z with homogeneous Dirichlet boundary conditions (in 
typical applications d = 1,2,3), where a = a{y) is a diffusion coefficient that depends on y, and 
where / G H~^{D) is fixed. The so-called affine case refers to a diffusion coefficient of the form 

a{y) = d + J2yj'^v (6) 

i>i 
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where a and the i'tpj)j>i are given functions from L°°{D). In this case, the yj typically range on 
hnite intervals which upon renormalization of a and ipj can be assumed to be [—1,1], Therefore the 
parameter domain is the infinite-dimensional box 

U ■.= [-!,If, (7) 

and well-posedness of the problem in V = HfD) is ensured for all y £ U by the so-called uniform 
ellipticity assumption 

^ |V’j(x)| < a(x) - r, xeD, (8) 

i>i 

for some hxed r > 0. For this model problem, convergence results have been obtained for polynomial 
approximations where = n, constructed by best n-term truncation of inhnite polynomial 

expansions, either of Taylor or Legendre type, that is, by retaining the n coefficients of largest norms 
in such expansions. 

A striking result, hrst established in [8] under the uniform ellipticity assumption, states that 
whenever the sequence (||'i/’j||L°°)j>i is summable for some 0 < p < 1, then such n term poly¬ 
nomial approximations converge with rate n~^ in L°°[U,V) where s := ^ — 1. For the Legendre 
approximations, an improved rate ri“®, with s ^ — i, is achieved in Lf{U,V,p) where y is 
the multivariate uniform measure, that is, in the mean-square sense if the yj are i.i.d. uniformly 
distributed variables. These results have been extended to a large range of linear or nonlinear para¬ 
metric PDFs [7, 9] where y is again ranging in the inhnite dimensional box U. They heavily rely 
on the holomorphy of the solution map y i—^ u{y) in each variable yj. 

1.2 Elliptic PDEs with lognormal coefficients 

In the present paper, we focus our attention on the so-called lognormal case for the elliptic PDF 
(5). In this case, the diffusion coeffient a is of the form 

a = exp(6), (9) 

where 6 is a random function of the form 

b = b{y) = '^yj'ipj, (10) 

i>i 

dehned from a given sequence {fjj)j>i of functions in L°°{D), with y = {yj)j>i a sequence of i.i.d. 
AA(0,1) variables. Thus, the parameter vector y now ranges over the unbounded domain 

U:=M.^. (11) 

We work with the usual product measure space given by 

= ( 12 ) 
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where B{U) = denotes the S-algebra generated by the Borel cylinders and 7 the tensorized 

Gaussian probability measure. 

We discuss further several conditions on the family {'4’j)j>i which ensure that the series in (10) 
converges almost surely in L°°[D). Under such conditions, b[y) is a Gaussian random variable 
with values in L°°{D). By the Lax-Milgram lemma, the associated weak solution u{y) of (5) with 
a = a{y) = exp( 6 (y)) is a random variable with values in U = Hq{D). 

One frequently used approach that leads to this framework is by starting with b = (b{x))x£D 
dehned as a centered Gaussian process over the domain D, with prescribed covariance function 

Cb{x, x') := K(b{x)b{x')^, x,x'&D. (13) 

One then obtains a representation of the form (10) by considering the Karhunen-Loeve decomposi¬ 
tion 

b = b{0 = (14) 

i>i 

Here are the L^-orthonormal basis of eigenfunctions of the integral operator with kernel Cb 

and are independent centered Gaussian variables. One then sets 

Vj ■■= and -01 := ■= 1E(^|). (15) 

Here again, conditions on {ipj)j>i are needed in order to ensure that the series in (10) converges 
almost surely in L°°{D). 

In this paper, we prefer to start with b dehned through (10) for a more general family ('ijjj)j>i. In 
particular we want to consider cases where the family (V’j)j>i in ( 10 ) is nonorthogonal and therefore 
differs from a normalized Karhunen-Loeve basis of b. One typical example is the case where b is 
a one-dimensional Brownian motion, which has a natural expansion in terms of the nonorthogonal 
Schauder basis. As shown later in this paper, this representation appears to be more efficient than 
the Karhunen-Loeve expansion in terms of the resulting best n-term polynomial approximation 
rates. 

In contrast to the affine model, the lognormal model is intrinsically stochastic: it reflects the 
situation where the diffusion coefficient a is allowed to become arbitrarily small or large (unlike in the 
affine case with the uniform ellipticity assumption), however with the probability being controlled 
by the Gaussian distribution. In this case, the relevant polynomial expansion is the orthonormal 
Hermite series, 

u='^UuH^, u^:= u{y)Hy{y)d'y{y), Hy{y) = W (16) 

■’U j>i 

where the univariate Hermite polynomials are normalized in L^(]R, dg), with g denoting the stan¬ 
dard Gaussian density. Note that is an orthonormal basis for the Hilbert space L‘^{U, M, 7 ). 

The coefficients Ui, are well dehned in V and the above series converges in the Bochner space 
L^([/, U, 7 ) associated to L^([/, M, 7 ) whenever u belongs to L^([/, U, 7 ). 


4 


Note that a function u \ y ^ u{y) in L‘^{U,V,'y) may also formally be viewed as a function 
{x,y) eT” u{x,y) := u{y){x) of both variables x G D and y and the inner product between two 
functions u,v ^ Lp‘{U,V,'y) has the expression 


{u,v) 



U JD 


Vrt(x, y) ■ Vv{x, y) dx d'y{y), 


(17) 


where the operator V is always meant in the x variable. In order to simplify notation, we shall 
systematically avoid mentioning the x variable, except when necessary, by simply writing the above 
integral as 

{u,v) = [ [ Vu{y) ■Vv{y)d'y{y). (18) 

JU JD 

The lognormal model and its approximation have been studied in various papers, e.g. [1, 5, 11, 
12, 16, 13, 17, 19, 21]. Sufficient conditions for the finiteness of moments E(||rt(y)||(>) at all orders 
0 < A; < oo have been established, either by assuming smoothness properties of the covariance 
kernel Cb or by assuming a summability property for the sequence (||'0j ||ioo(£)))j>i, see for example 
[5, 11, 19]. Since 

®"(ll^(y)lly) ~ (19) 




the finiteness of the second moment is a necessary and sufficient condition for the summability 
of the sequence (||rtjy||y),ygj- and the convergence of the Hermite series in L^(t/, 17,q). However, in 
order to prove convergence rates for the best n-term truncation 


rtAn •“ ^ ^ UyHy, (20) 

of (16), where A„ denotes the set of indices corresponding to the n largest ||ttjy||y, we need to study 
the summability properties of (||tii/||y),ygj- for p < 2. For example, if this sequence is proven to 
be iP summable, then a direct application of Stechkin’s lemma [10, 8, 7], implies 

Ik “ ^AnllL 2 (f/,y,^) = ^ <C{n + l) (21) 

where s := | - 1 and (7 := \\i\\uu\\v)u<^D\\^p(^jry 

The first and currently the only available result concerning summability of Hermite coeffi¬ 
cients, in the infinite dimensional framework, has been established in [19] and reads as follows. 


Theorem 1.1 For any 0 < p < 1, if (j||kil|i/°°)j>i ^ ^^(N) then {\\uy\\v)v^T £ ^^(7^)- 

The Hermite coefficients studied in [19] are actually different from those defined in (16), since 
the authors are interested in a specific approximation process, namely the Galerkin method. One 
specihc difficulty in this approach is the fact that the bilinear form 

B{u,v)=[ [ a{y)Vu{y) ■Vv{y)d'y{y), (22) 

JU JD 
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is not continuous and coercive on L'^{U,V,^). This particular issue has been dealt with in papers 
such as [5, 13, 16, 12]. The approach from [16], which is also used in [19], is based on a modified 
bilinear form where 7 is replaced by an auxiliary Gaussian measure. The Galerkin error can then 
be estimated by the tails of modified Hermite coefficients corresponding to a further auxiliary 
Gaussian measure. However, inspection of the proof of Theorem 1.1 in [19] shows that it also holds 
for the standard Hermite coefficients. Note that this issue is specific to the Galerkin method, and 
does not occur in pseudo-spectral or least-squares methods. 

In the present paper, we only consider the classical Hermite coefficients in (16), since our main 
interest is to obtain general polynomial approximation results for the solution map, that may further 
serve as a benchmark for any numerical method based on such approximations. 

The above Theorem 1.1 is comparable to those obtained in [ 8 ] for Taylor and Legendre coefficients 
in the affine case, except for the appearance of the additional factor j in front of ||'0j||L°°) which 
makes the summability assumption more restrictive. On the one hand, it is quite remarkable that 
such a result exists, since the complex variable arguments used in the affine case, based solely on 
holomorphy of the solution map, are not sufficient to control the Hermite coefficients. The proof 
of the above result in [19] relies, instead, on estimates for the mixed partial derivatives of u up to 
some finite order m related to p. On the other hand, the type of assumptions used in this result 
indicates that it is not in the sharpest possible form as we now explain. 

First, the above result does not cover summability for 1 < p < 2, which one might expect to be 
obtainable under weaker assumptions. Second, and more importantly, the summability conditions 
imposed on the ||x,c 3 o in these results become quite strong and artificial in the case when the 
supports of these functions do not overlap too much. 

As a relevant example for this second point, consider the case where the ^ire a wavelet 

basis on the domain D. In this case it is more natural to index such bases according to (V’A)AeV! 
where A is a scale-space index following the usual terminology, such as in [ 6 ]. With the notation 
/ = |A| for the scale level, there are wavelets at this level and each of them has support of 

diameter 0{2~^). The supports of wavelets at a given scale I have finite overlap, in the sense that 
any x £ D is contained in the support of at most M wavelets of level I where M is independent of 
X and 1. We assume that the norms of the wavelets only depend on the scale level, that is, 

||^a||l°° = Q) 1^1=^- (23) 

It is well known that the geometric rate of decay of the wavelet contributions, as the scale level 
grows, reflects the amount of smoothness in the expansion (or the smoothness of the correlation 
function in the case of a random series). It is thus natural to study the situation where q is of the 
form 

Cl := C2-^\ (24) 

for some given a > 0. Then, it can be checked, for example using the arguments from [5], which 
are recalled further in §2, that for arbitrary a > 0 and C > 0, the bound (24) implies that the 
solution map has bounded moments E(||n(y)||y) for all 0 < A: < 00 . Indeed, this follows from 
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the Lax-Milgram lemma and the fact that E(exp(A:|| 6 (y)||i;,oo)) is finite. The particular case k = 2 
implies that the Hermite coefficient sequence (||rtjy||y),ygj- belongs to provided only that 

a > 0. However, if we want to use the above mentioned result from [19] to prove summability of 
this sequence for values of p < 2 , the appearance of the factor j, suddenly, requires the very strong 
constraint a > 2d which implies a strong smoothness condition on the diffusion coefficients. 

The above wavelet example reveals a gap in the currently available analysis: summability can 

be obtained under mild assumptions on the smoothness of the diffusion coefficient, while proving 
summability for p < 2 by the existing results immediatley imposes much higher smoothness (in 
the sense of the required decay of as j —)• oo). The goal of the present paper is to propose 

a sharper analysis, which removes this gap. Let us mention that the same gap occurs, in a slightly 
less pronounced form, for the Legendre coefficients in the affine case, an issue which we recently 
adressed in [ 2 ]. 

1.3 Main results and outline of the paper 

The main result of our paper for the model elliptic equation (5) with lognormal diffusion coefficient 
is the following. 

Theorem 1.2 Let p < 2 and let q = q{jp) := Assume that there exists a positive sequence 

(Pj) such that 

snp'^Pj\'iljj{x)\ < oo. (25) 

and 

G ^^(N). (26) 

Then, the solution map y i— u{y) belongs to L^{U, V, 7 ) for all 0 < k < 00 . Moreover, (||tti/||y)jygjr g 
In particular, best n-term Hermite approximations converge in L^(t/, H, 7 ) with rate n~^ 
where s = J - i = i. 

This theorem gives a significant improvement over Theorem 1.1 in the case where the functions 
ijjj do not overlap too much. First of all, in the case of disjoint supports, the first condition in 
Theorem 1.2 is met when pj^ := Therefore, it implies that, for all 0 < p < 2 the 

H summability of (||rii/||y)i/gj- follows from the assumption that (HV'j l|L°°)j>i is summable for 
q = q{p) := Note that q{p) > p for any p > 0 and that 

lim q{p) = + 00 , (27) 

p —^2 

which shows that almost no decay of (||V'j||i,°°)j>i is required as p gets closer to 2 . Secondly, as is 
shown later in this paper, we can also use Theorem 1.2 to treat the above mentioned wavelet case, 
and obtain summability results with p < 2 for any smoothness index a > 0. Finally, even in the 
case of arbitrarily supported xfj, we establish as a corollary of Theorem 1.2 that the H summability 
of (||Ui/||y),ygj- follows from the summability of {\\'4’j\\L°°)j>i with q := which again represents 
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a significant improvement over the condition in Theorem 1.1, since this summability follows by 
Holder’s inequality from the summability of {j^\\'4’j\\L°°)j>i for some /3 > ^. 

Let us mention that a result analogous to Theorem 1.2 has been established in [2] for Taylor 
and Legendre coefficients in the affine case, under the additional uniform ellipticity assumption. 
A common feature with the present paper is that one key ingredient of both proofs consists in 
establishing certain weighted estimates for the coefficients ||njy||y which further translate into P 
estimates by Holder’s inequality. However, establishing such weighted estimates for the Hermite 
coefficients in the present lognormal cases uses completely different techniques than those in [2]. In 
particular, similar to [19], we only rely on the study of mixed partial derivatives d^u for a limited 
order ||/i||£oo < r, while the estimates in [2] exploit all orders ^ J-. Let us also mention that, while 
we focus on the diffusion equation, inspection of our proofs reveals that our main results can be 
extended to other types of linear elliptic or parabolic PDEs with similar lognormal dependence of 
the coefficients. 

The rest of our paper is organized as follows. We begin in §2 by revisiting conditions which 
ensure that y u{y) is a measurable map with values in V and with finite moments E(||tt(y)||y) 
for all 0 < A: < oo. We introduce a sufficient assumption in terms of the convergence in L°°{D) of 
the series Pj\'4’j\ for a positive sequence {pj)j>i which satisfies the summability property 

^exp(-p2) < oo. (28) 

i>i 

This is a very weak condition on the sequence {pj)j>i and in particular always holds under the 
assumptions in Theorem 1.2. 

As a first step in the proof of Theorem 1.2, we relate in §3 the norms ||u,y||y of Hermite coefficients 
with ||i 9 ^u||^ 2 ([/ y,|,) for relevant values of y,. This leads us to an identity between certain weighted 
norms of both quantities. As mentioned above, while all y G J-" are considered for the Hermite 
coefficients, only limited order ||/i||£oo < r are considered for the partial derivatives. As a second 
step, we obtain in §4 bounds on the previously introduced weighted norms of the \\d^u\\]^ 2 (^ijy^^-^ 
under the assumption that Pj IV’i(®)l bounded by a relevant constant. Finally, we combine 

these ingredients in §5 to complete the proof of Theorem 1.2. 

In § 6 , we derive various consequences of Theorem 1.2 corresponding to the different cases outlined 
above for the functions (i) disjoint or finitely overlapping supports, (ii) wavelets, (iii) 

arbitrary supports. We discuss, in all three cases, the improvements over Theorem 1.1. 

We conclude in §7 with an interesting observation which follows from our analysis, in the par¬ 
ticular case where 6 is a Brownian motion. For this case we compare the convergence of the best 
n-term truncation of Hermite series when using either the Karhunen-Loeve representation or the 
Schauder basis representation for b. The first representation satisfies by construction an orthog¬ 
onality between the ijjj, while the second one does not. Due to the limited amount of smoothness of 
b, the assumptions of Theorem 1.2 do not hold for the Karhunen-Loeve representation, so that no 
algebraic convergence rate can be established for the resulting polynomial approximation with our 
currently available techniques. In contrast, we can obtain an algebraic convergence rate when using 
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the Schauder representation, by exploiting the local support properties of the basis functions. This 
hints that, in relevant practical cases, the Karhunen-Loeve representation might be suboptimal in 
terms of the resulting polynomial approximation rates for the solution map y i—)• u{y) in L?‘{U, 1 ^, 7 ). 

Let us stress that our results only quantify the approximability of the solution map. They should 
therefore be viewed as benchmark for concrete numerical methods. The development of numerical 
methods that provably meet such benchmarks will be the object of further investigation. 

2 Measurability and integrability of the solution map 

The hrst question we investigate is what conditions on {'4’j)j>i guarantee that the solution map 
y I—>■ u{y) is in L^{U, V,7), i.e. it is measurable with values in V and satishes 

^{\Hy)\\v) = [ \\u{y)\\^ d-yiy) < 00 . (29) 

Ju 

We recall that a function from a measurable space to a Banach space B is measurable (also sometimes 
called /r-measurable, or Bochner-measurable, or strongly measurable), if and only if it is the almost 
everywhere pointwise limit of a sequence of simple functions. 

Recall that u{y) is the weak solution to the diffusion equation 

-div(a(y)Vu(y)) = /, (30) 

on the bounded domain D with homogeneous Dirichlet boundary condition, for a fixed / £ 
were 

a{y) := exp(6(y)). (31) 

The solution u{y) is well dehned as an element oiV = Hq{D) provided that Vjipj defines a 

function b{y) G In such a case, by the Lax-Milgram lemma, we have 

\\u{y)\\v <C\\a{y)~^\\L-o <Cexp{\\b{y)\\L^), C := ||/||^-i. (32) 

This motivates the study of the finiteness of exponential moments of || 6 (y)|| 2 ,°o which we formulate 
for any given 0 < fc < 00 as follows. 

Property M^: The map y 1 —)■ b{y) is measurable from U to L°°{D), and the exponential mo¬ 
ment E(exp(/c|| 6 (y)|| )) is finite. 

According to (32), we hnd that Property implies the finiteness of the A;-th moment E(||ri(y)||y). 
Several types of conditions on b have been introduced in the literature [5, 11, 19] which are sufficient 
to guarantee that Property Mfc holds. 

One first approach requires that b = {b{x))x£D is a Gaussian random held with some minimal 
smoothness, in the sense of a smoothness assumptions on the covariance kernel Cf,. Indeed, if Cb 
belongs to the Holder space C^{D x D) for some /3 > 0, that is, there exists C > 0 such that 

\Cb{x^x) — Cb{z,z')\<C{\x — z\^+ \x'— z'\^)^ x,z,x',z'eD, (33) 
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then, following the argument of Proposition 2.1 in [5], we first find by using the Kolmogorov conti¬ 
nuity theorem that there exists a version of b with trajectories almost surely in C°^{D) for a < (3/2. 
By application of Fernique’s theorem, see Proposition 2.3 of [5], we then find that 

E(exp(A||6||icx,)) < E(exp(A||6||^„)) < oo, (34) 


for some sufficiently small A > 0. This implies that E(exp(A;||6||Loo)) is finite for all 0 < /c < oo. 
Note that when h has the representation (10), its covariance kernel is then given by 

Cb{x,x') = '^iljj{x)'4)j{x'), (35) 

i>i 

and therefore we may analyze the smoothness of Cf, through that of the individual ipj. A strategy 
to establish that b with representation (10) is almost surely in a Holder space for some k > 0 
is discussed in [18, 11]. It assumes that the ijjj are individually in C" for some a > 0, and that for 
some 0 < (5 < 2, 

W'fpjWl^ < oo and ^ < oo, (36) 

i>i i>i 

where Cj is a constant such that |'0j(®) ~ ^ Cj\x — x'|“. Then, using the Kolmogorov 

continuity theorem, one establishes that b belongs almost surely to for k < a5j2, see Corollary 
7.22 in [11] (one can also easily check that (36) implies (33) with (5 = 5). 

A second approach for guaranteeing Property does not assume any smoothness on the ipj, 
but instead the summability property 


^ IIV'jIIl'- < oo. (37) 

i>i 

Then, as observed in [19], we find that for all 0 < /c < oo, 

E(exp(A:||6(y)||Loo)) <E(exp(^k'^ |yjHIV’illL°°)) = 

i>i l>i 

since for a standard Gaussian variable t and a positive quantity s one has 

2 r-l-oo / ^2 \ 2 P / §2 2 A 

E(exp(»|t|)) = ^y__ “»('*- 2 )"* = '^757-. 

Neither of these two types of conditions imply each other. Indeed, on the one hand, if {'4’j)j>o 
is a smooth wavelet basis of L^{D) ordered from coarser to finer scale, smoothness of Cb holds 
whenever ||'i/’j||L°° < Cj~^, which is not sufficient to ensure that (||V’il|L°°)i>i belongs to £^(N) if 
^ < 1. This shows that (33) may hold while (37) may fail. Also note that, for the same reason, 
the criterion (36) is sufficient but not necessary for (33) to hold. On the other hand, if the ipj are 
discontinuous, we cannot hope that Cb has Holder smoothness while (||^j ||i,°°)j>i could belong to 
.£^(N), which shows that (37) may hold while (33) and (36) may fail. Note that the exclusion of 
discontinuous V'j is problematic when one needs to model sharp interfaces in the diffusion media. 
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We also remark that conditions (36) and (37) are artificially strong when the supports of the ipj 
do not overlap, or only partially overlap. This motivates us to introduce a third condition which 
better takes into account the support properties of the i/jj , and does not enforce them to be Holder 
continuous. 

Assumption A: There exists a strictly positive sequence p = such that the series PjIV’il 

converges in L°°{D) and 

^exp(-p|) < oo. (40) 

i>i 

Note that since the pj are strictly positive, the summability of the exp(—p^) also implies that 

Po ■= inf Pj > 0. (41) 

We shall make use of the following elementary lemma. 

Lemma 2.1 If {aj)j>i is a sequence of numbers from [0,1[ that belongs to -^^(N), then, there is a 
constant c > 0 and Mq > 1 such that 

OO 

M>Mo. (42) 

i=i 

Proof: We first consider the case when Oj < 1- In this case, (42) follows with Mq = 1 because 
||(Q:j)j>i||£M < ||(aj)j>i||£i, for any M > 1. Hence, with c= — In ||(aj)j>i||£i, we arrive at (42). To 
prove the general case, we first choose K so that — I/^- ff ^0 > 1 is sufficiently large 

then 

OO K 

E “f“ = E “f“ + E “f“ S 3/4. (43) 

j=l j=l j>K 

where we used the case already proven for estimating the second term. Since 
know from the comparison of norms 

OO 

< (3/4)^/^°, M > Mq. (44) 

1=1 

and the result follows. □ 

We now prove that Assumption A guarantees the measurability and finiteness of the exponential 
moments of the map y e-)• b{y). 

Theorem 2.1 Assumption A implies that Property holds for any 0 < k < oo. 
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Proof: Let us consider a sequence {pj)j>i such that (40) holds. For any t >0, the complement 
of the event 

£t := {y : suppJ^lyjl < t}, 


i>i 


has measure 


liSf) ■ \yj\ > tpj} < 

i=i ” j=i 




where we have used the univariate Gaussian tail bound 


/ 

'|s|>B By^TT 


e 2 B > 0. 


Therefore, application of Lemma 2.1 with aj = exp(—/9|) shows that 

7(^f) < t>to, 

for certain constants C, c, to. In particular the event 

£ := Ut>o£t = {y : supp~^\yj\ < oo}, 


i>i 


has full measure, i.e. 


j{£) = 1 . 

We next introduce for J > 1 the truncation 

J 

bj{y) = '^yj'>Pj- 

i=i 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


The mapping y i—>• bj{y) is measurable from U to L°°{D) since it is a continuous L°°-valued function 
of the variables {yi ,..., yj). We observe that if Assumption A holds, then, for any y € £, we may 
define b{y) as the limit in L°°{D) of bj(y) since 


\\bj{y)-b{y)\\L^ < ^ Pjl'i/'jl suppj^\yj\^0 as J —)■+ 00 . 

j>J 


(52) 


Thus bj(y) converges to b{y) in L°°{D) for almost every y £ U. Therefore y e-?■ b{y) is also a 
measurable mapping, that is, b{y) is a random variable with values in L°°{D). 

In addition, for any s > 0 and y € we may write 


Hence, for t> 1, we have 


Lcc < Cas, Ca ■■= II jOj#j 
i>i 


P{t) := 7{y : ||6(y)||L- > t} < j 


(53) 




(54) 
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where we have used (46). From Assumption A, we know that the last sum in (54) is finite when 


/ ^2 p2 \ 

t = to := V2Ca- Hence, applying Lemma 2.1 with aj = exp(— 

and using the fact that 

P{t) < 1 for all t > 0, we find that 


P{t) < Ce-^^\ 

(55) 

for suitable c, C > 0 and for all t > 0. Therefore 


^+00 


E(exp(fc 6 (y) L°°)) = / kex.p{kt)P{t)dt < 00 , 

Jo 

(56) 

for all 0 < A: < 00 . 

□ 


If we next define uj{y) as the weak solution of (5) with a = aj{y) = exp( 6 j(y)), where bj is as 
in (51), we find that the mapping y i—)• uj{y) is measurable from U to V since it is a continuous 
f/-valued function of the variables (yi, ... ,yj). This continuity stems from the following classical 
stability estimate: two weak solutions u and u of (5) with diffusion coefficients a and d, respectively, 
satisfy 

\\u-u\\v <C\\a-d\\L^, C:= . J*^**^! -(57) 

This estimate also shows that the convergence of bj{y) towards b{y) in L°°{D) for each y & £ 
implies the convergence of uj{y) towards u{y) in V for each y £ £, which shows that the mapping 
y e-)• u{y) is measurable as an almost everywhere limit of measurable mappings, that is, u{y) is a 
random variable with values in V. By application of (32), we thus obtain the following result. 

Corollary 2.1 Assumption A implies that y i— u{y) is measurable with values inV = Hq{D) and 
that E(||rt(y)||y) is finite for any 0 < k < oo. 


Remark 2.1 Assumption A is almost necessary for Property to hold in the case where the 
supports of fij do not overlap, in the sense that Property then implies that is 

uniformly bounded. In this ease, the uniform boundedness o/PilV’il for o, sequence {pj)j>i 
which satisfies ( 40 ) can be equivalently expressed as follows: there exists Cb > 0 such that with 


bj := 11 11 /^ 00 , one has 


^exp(-^ 


j>i 


< oo. 


(58) 


Indeed, since Ca '■= = supj>;^pj5j, we have (58) with Cb = . We then have 

11 L'^ 


i{y ■ \\b{y)\\L-°{D) < =n ^{2/ • - n 


i>i 

for any t > 0, where we have used the lower bound 


f>i 


2B 


1 - 


2t 


bj\/^ 


exp 


/ P 


'ld>B 


g{s)ds> / g{s)ds> _ exp(— 

Jb<\s\<2B v27r 


(59) 


(60) 


It follows that if (58) does not hold for any Cb > 0, then 'y{y : ||ft(y)||L°°(D) < = 0 for all t > 0, 

which means that \\b{y)\\B^{D) = 00 with probability 1 and therefore Property does not hold. 
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Remark 2.2 Assumption A always follows from the assumptions of Theorem 1.2. Indeed, if 
ipj^)j>i G for some q < oo, one can find a positive sequence such that linij^+oo dj = 0 

and such that the sequence defined as loj := 6jPj satisfies (exp(—G ^^(N). We then obtain 
the validity of Assumption A with the sequence (ujj)j>o since 


j>J 




—)• 0 , 


as J —)• + 00 . 


(61) 


3 Hermite expansions and partial derivatives 


We are interested in the summability properties of the multivariate Hermite expansion (16). By 
Parseval’s identity, the integrability of u{y) with respect to the Gaussian measure 7 implies the 
summability property 

\\uu\\y < 00 . (62) 

In order to obtain polynomial approximation results in the mean-square sense, or equivalently, in 
L^(M^, 1 /, 7 ), we need to establish the summability 


\\uu\\v < OO’ (63) 

v&T 


for values of p < 2. For this purpose, we establish L^-integrability properties of the partial deriva¬ 
tives of u. For any H-valued function y i— wiy) defined on and f G N, the derivative dy.w is 
defined as the limit, as h —)• 0 , of the difference quotient 


wiy T hcj) - w{y) 
h 


(64) 


provided this limit exists in V, where 

Cj .— ((^i,j)i>l (65) 

is the Kronecker sequence of index j. Higher derivatives d'^w are defined inductively. For v d IF, 
we use the notation 

d'^w = dy^j^w- ( 66 ) 

j>i 

When, Uj = 0, the operator dy^ is the identity operator and so the above product has only a finite 
number of factors which are not the identity. 

We also use the standard notation supp(z^) := {j : vj 7 ^ 0}, as well as 

:= II^^IIh = (67) 

i>i 


and 


z/! := z/j!, z/ G F, 

1>1 


( 68 ) 
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with the convention that 0! = 1. Furthermore, we use the combinatorial notation, 


n 

i>i 




with the convention that 


:= 0 , if m > n. 

In analogy with y = {yj)j>i, we set 'ip = i'4’j)j>i and use the notation 

i>i 


( 69 ) 


(70) 


(71) 


Finally, with < denoting the componentwise partial order between multi-indices, we define 

■.= {v £ T ■. V < [I and v / y}, y £ F, (72) 

which is always a finite set. The following lemma proves that the derivatives d'^u{y) are elements 
of V and gives a recursive way of computing them. While it is implicitly contained in the proof of 
Theorem 3.1 in [19], we give it here in an explicit form for completeness. 

Lemma 3.1 For any y such that ||6(y)||L°° < oo, and any nonzero v £ F, the partial derivative 
d'^u{y) exists and is the solution of the variational problem 

[ a(y)Vd>^u(y) • Vu = - V f a{y)Vuiy) -Vv, v£ V. (73) 

Jd Jd 

Proof: The proof of this lemma is, in principle, the same as the proof of Theorem 4.2 in [ 8 ]. The 
only difference is that in the theorem of [ 8 ] the assumption is that a is affine in y, whereas now we 
assume a = with b affine in y. In going further, we indicate only the necessary changes caused 
by the different diffusion coefficient. 

One begins by establishing the validity of the theorem when y = ej, for any j £ N, by considering 
the functions ph{y) ■= ^ for \h\ < 1. Following the reasoning in [ 8 ], one derives that p/i 

satisfies 

a{y)yph{y) ■ Vu = Lh{v), for all v £V, 


L 


where Lh is the linear functional on V given by 

g/ifo _ 

Id h 


r _ I 

Lh{v) ■■= - - 7 - a{y)Vu{y + hej) ■ Vv. 

Jd ^ 


Using (57), one next proves that L/j —)• L in V' where L is the linear functional 

L{v) ■■= - a{y)'ipjVu{y) ■ Vv. 

Jd 


(74) 


(75) 


(76) 


15 




Therefore, taking a limit as /i —)• 0 in (74) proves the theorem for the case /r = Cj. The case of any 
general /r G now follows by a recursive application of this result. □ 

The recursive formula (73) will be crucial for obtaining sharp estimates for the partial derivatives 
d^u{y). For the time being, we use it to establish the following crude bound 

||9'''(^(y)||l/< C'exp((2|/r| + l)||6(y)||Loo), /r G J’, (77) 

where the constant depends on y, on the ipj, and on /. In view of Property M^, this guarantees 
the finiteness of all moments E(||cl^tt(y)||y) for all 0 < A: < oo. 

For the proof of (77), for each y such that ||&(?/)||l°o < oo, we introduce a y-dependent equivalent 
norm on V, 

\Ml(y) ■■= o(y)|Vup dx, v£V. (78) 

For such y, we have the norm equivalence 

exp(-||6(y)||Loc)||u||^ < lblla(j/) ^ exp(||6(y)||Loc) ||r;||^. (79) 

Using this in conjunction with (73) and v = d^u{y), we obtain 

<^ (U( ||%|i;o‘'')||a'u(!/)lUteil|a'‘u(!/)lUtei 

I'&Sij, jGsupp{y,—i/) 

<C^e2 lW^)ll^-||a^u(y)||y \\d''u{y)\\v, 

with a constant Cy > 0. Applying this recursively, starting from the estimate (32) for ||tt(y)||r/, we 
obtain (77). 

The following result relates the norms of the mixed derivatives of u up to some given order 
r with weighted norms of the Hermite coefficients. 

Theorem 3.1 Let r > 0 be an integer and let p = he a sequence of positive numbers. 

Assume that Property hold for all 0 < k < oo. Then 

^ = '^K\K\\v, ( 80 ) 

\\uh^<r 

where 

E (81) 

The weights by can also he expressed as 

i'-=n(E(7)d')^ (82) 

j>l 1=0 ^ ' 


16 


Proof: With the normalization in L‘^(K,dg) that we have chosen for the univariate Hermite poly¬ 
nomials, we have the Rodrigues formula 


Hn{t) = 


(_l)n^(n)(i) 


(83) 


a/uJ git) 

As a consequence, for m < n and sufficiently smooth univariate functions v with square-integrable 
derivatives, integration by parts yields 


Vn-= / vHngdt = 

Jm. 

By Parseval’s identity, we thus have 

^ [ \v^^^\‘^gdt = V 

Jr 


dt = t^ [vgi-Ut = \[^^ 
Jr V n! 


J^^Hn-mgdt 


ni 


n>m 


m!(n — m)\ 


Vn = 


E 

n>0 




(84) 


(85) 


where in the last step we have made use of our notational convention (70). 

For z^, // E with // < zz, applying the same reasoning to the tensor product Hermite coefficients 
Uiy, we obtain 

Uu = (U, Hy) = Hy-fj) . (86) 

Concerning the integration by parts, note that, as a consequence of (77), Property allows us to 
proceed as in (84). By Parseval’s identity, summation over zz E with v > g thus gives 


1 







b/ l|2 



(87) 


where we have again used (70). Multiplying this identity by and summing over g with ||/z||£°o < r 
gives 

^ [ \\d^^iy)\\vd-fiy)= = ( 88 ) 

that is, (80) with by as defined in (81). The equivalent expression of by, given in (82), immediately 
follows by factorization in (81). □ 


4 Estimates of the partial derivatives 

fn view of Theorem 3.1, estimating certain weighted norms of the sequence (||ui/|| v")jygj^ amounts 
to estimating related weighted norms of the sequence 

([ \\d^uiy)fy djiy)) , (89) 

\d U J ||/.i|hoo<r 

which we next address. 
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4.1 Estimates for fixed y 

Recalling the norm || • ||a(y) defined in (78), we first fix y and estimate the relevant weighted 
norms for the sequence (||d'^';^(y)||a(y))||^||^^<^- 

Theorem 4.1 Let r > 1 be an integer. Assume that there exists a positive sequence {pj)j>i such 
that 

sup '^Pj\i/jj{x)\ =: K < Cr := — 1 = ■ (90) 

x&D 

Then, there exists a constant C that depends on K and r, such that 

E ^ lia"'‘fe)lli(,)<CM!;)lli(,,, (91) 

llMl|oo<r 

for all y satisfying \\b{y)\\ < OO. 


Proof: For k >0, we define := {p £ IF : \p\ = k, ll/illr;,^ < r} and 


We prove that 


/iSAj. 




for a fixed (f < 1. Since cJo = ||u(y)||^^^^, the theorem follows from this by summing over k. 
We introduce the notation 

ip-uy. ’ 

where := OjM From (73), recalling the notation (72), we have 

p'^\Vd'^u{y)\ pk'\Vdk'u{y)\ 

-pp^ujuyy^ 

and thus, by Cauchy-Schwarz inequality. 


Ufc S 


^ [ Y1 ^il^F)a{y)- 
Jd .trf 


^ilk,u)a{yy- 




For p £ Ak, define -.= {u £ S^-. \p — v\ = A}. We have 

k 


= Y Y 


(92) 


(93) 


(94) 


(95) 


Y ^iLF)a{y)- -j- . (96) 

,,^o A- / 


(97) 




£=i veS,j,t 
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Now, from (90), we have 


e(/r, u) < //2 




< E 


(u, — uV. 

OO 0 1 

-I 


T\ 


= //2 


P 1 


i=i 


Inserting this into (97) we find that 

^ e(/i, v) < - 1 < -1 = 1 . 


( 98 ) 


Inserting this into (96), and applying the Cauchy-Schwarz inequality, we obtain 
CTk < e{n,u)a{y) 


\p''Vd''u{y)\‘^V^^ / , Jp/^V5^^r(y)|2^ 

a{y)- 






- / f E E 


2^ a{y)- 




fiGAfg uGSfj, //EA/j 

We treat the first factor in the last integral by interchanging summations in p and v. For this 
purpose, we introduce, for any i < k — 1 and ly € A^, the set 


Ry,k ■— ■ k' G < 5 ^}' 

By a similar argument as above for the set 5*^/, we find that 


E - 




1 


fJ'GRi^^k 


ik-i)l 


K 


k-£ 


(99) 

( 100 ) 


It follows that 


o-fc < 


'D 


k-l 

E- 

e=o 


.(k-i)l2_ 


E -iy) 


|p^V5^M(y)|^ 


2 \ 1/2 / 


{k-i)\ 


u&Ai 


V\ 


(E “(y) 

/iGAfc 


\pk-'S/ dku{y){" 


2 \ 1/2 


If we now apply the Cauchy-Schwarz inequality on the integral we obtain 

k-l 


1=0 ^ ’ 


1/2 


^k 


1/2 


In other words, 


k-l 


^k<Y 


1 




i=o 


{k-ey. 


( 101 ) 

( 102 ) 

(103) 


Now pick (5 < 1 such that K/6 < Cr- By induction, we prove that ak < for all k > 0. This is 
clearly true for k = 0 and assuming it has been proven for A: — 1, we find 


(^k < <70 


' E ^ ^ 0 ^" E ^ - 1 ) = lTo<5^ (104) 


{k-iy 


1=0 ' ^ 1=0 
which advances the induction. This establishes (93) and completes the proof of the theorem. □ 
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4.2 Integral estimates 

Combining (79) with Theorem 4.1, we obtain the following result. 


sup 

x&D 


Theorem 4.2 Let r > 1 he an integer. Assume that there exists a positive sequence {pj)j>i such 
that 

” Pj\^Pj{x)\ =:K<Cr := (105) 

■j>i 

Assume in addition that Property M 4 holds. Then we have 

^ [ \\d^u{y)\\‘^ d-f{y) < 00 . (106) 


Proof: We apply (79) with v = d^u{y) and obtain 

||5^M(y)||^ < exp(||6(y)||Loc)||0^u(?/)||2(^). 

This gives 


(107) 




^ [ \\9^u{y)\\l dj{y) < j exp(|| 6 (y)||Loc) ^ 

11.00<r 


D <r 


< 


< 


< 


C'i^exp(|| 6 (y)||Loc)||u(y)||^(j^) d-f{y) 

Cl [ exp( 2 || 6 (y)||i,oo)||u( 2 /)||^(i 7 (y) 

Ju 

C 2 [ exp(4||6(y)||i,oo) (i 7 (y), 

Ju 


(108) 


where the second inequality uses Theorem 4.1, the third inequality uses (79) again, and the last 
inequality uses (32). From Property M 4 , the last integral is finite and the proof of the theorem is 
complete. □ 


5 Summability of Hermite coefficients 

According to Theorem 3.1, the conclusion of Theorem 4.2 also gives the weighted summability 
estimate 

^ ^ 1 / 11 ^* 1 /||y < 00 . (109) 

Using Holder’s inquality, for 0 < p < 2 and q := we thus have 

E ii“^iil s (E oAuAlf'iY: i-Af (110) 

uGT vGT 

The following result allows us to control the second factor. 
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Lemma 5.1 Let 0 < p < 2 and q := and let {pj)j>i be a positive sequence such that 


(Ill) 


Then, for any positive integer r such that < p, the sequence [by) defined in (81) for this r 
satisfies 

< oo. (112) 


E 


Proof: The sum under consideration can be rewritten in factorized form as 

= En(E(?)d')'”'' = nE(E 

u&T j>l 1=0 ^ ^ 1=0 

provided we can show that the product on the right side is finite. Now, 


^ ' 21 
l,Pj 


(113) 


E(E 

n>0 1=0 


I 1^^ 
with the constant 


-Q /2 

PT} < E 

n>0 ■- 


\ 2(nA'r) 
n Ar)^^ 


-1-9/2 


< I + Pj'^ + ■ ■ ■ + Pj^ + Cr,qPj (114) 


^'■-9 •“ E 


n>r 


-9/2 


(r!)*^/^ E [(^ + 1) ■ ■ ■ (^ + ’")] 


-9/2 


(115) 


n>0 


We have Cr,q < oo if and only if q > 2/r, which holds precisely under our assumption p > 

Since pj oo, there exists J such that pj > 1 for all j > J. For such j, we can bound the right 
side of (114) by 1 + {Cr,q + r — l)p~'^. Hence, returning to (113), for a finite constant C depending 
on w, we have 


E + {Cr,q + r - l)pj'^) , 

ly&T j>J 

where the product on the right side converges because X]j>i pj'^ < oo. 


(116) 

□ 


We are now ready to prove our main result. 


Proof of Theorem 1.2: By our assumptions, we are given a positive sequence {pj)j>i such that 

2p 


sup E Pj\fij{x)\<oo and {pj ^)j>i e q 


xGD 


i>i 


2-p 


(117) 


We choose r as the minimal integer satisfying < p. Since the sequence pj can be rescaled 
by multiplying by an arbitrary constant without affecting (117), we can assume without loss of 
generality that 

„ In 9 

(118) 


sup E/’j#j(a^)l < 


x£D 


i>i 
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Since, as noted in Remark 2.2, {p- ^)j>i G ^'^(N) implies Assumption A, by Theorem 2.1 we conse¬ 
quently also have Property for all k. We can thus combine Theorems 3.1 and 4.2 to obtain 

'^K\\uu\\v= y] ^ / Wu{y)\\vdi{y) <oo. (119) 

v&F ||At|k°o<r 

Using this together with Lemma 5.1 in (110), we arrive at the conclusion. □ 

Remark 5.1 As stated in Theorem 1.2, application of Stechkin’s lemma gives that the best n-term 
approximation polynomials 

UAn ■= y] UuH,,, (120) 

obtained by retaining the indices of the n largest satisfy the estimate 

IW - UAj\L2(^uy^^) < n“^, (121) 

where s := ^ ~ \ There is, however, a more direct and constructive way of retrieving this 

convergence rate, namely taking instead A„ to be the set of indices corresponding to the n smallest 
values of the weights by which appear in (109). We then directly obtain that 

1 i‘2i 

\\u - UAjL2(U,Vn) < sup bf^^‘^('^by\\uy\\l) < <+ 1 , (122) 

where {dl^)n>i is the decreasing rearrangement of the sequence {by ' )ueT- As seen in Lemma 5.1, 
this sequence belongs to 1^{T) which implies that d* < n~^ with s := K 


6 Examples 

In this section, we present several examples of applications of Theorem 1.2 corresponding to different 
support properties of the {fjj)j>i- In each case, we discuss which range of summability of the 
sequence (||V’j lli>i)i>i implies summability of the sequence (||rii/|| v")jygj^ for some p < 2. 

6.1 Finitely overlapping supports 

We say that the family {'ipj)j>i has finitely overlapping supports if and only if there exists an integer 
AL such that for every x £ D, 

#{j : V',(x)/0}<M. (123) 

One example with M = 1, which corresponds to disjoint supports, is the set of characteristic 
functions 

'i(j = CjXDj, (124) 

with some normalizing factor Cj, when {Dj)j>i is a partition of D. Another example with M > 1 is 
the set of Lagrange finite element basis functions of a given order k > 1, associated to a conforming 
simplicial partition of D. 
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For such families, we find that the choice 


Pi ^ — 




(125) 


yields 


sup E Pj\ipj{x)\ < M, 


xGD 


(126) 


i>i 


and therefore condition (25) in Theorem 1.2 is satisfied. We thus obtain the following immediate 
corollary. 

Corollary 6.1 Let {'ilij)j>i be a family with finitely overlapping supports, and let 0 < p < 2. If 
{\\'4’j\\L°°)j>i belongs to for q = q{p) := then (||rtjy||y)i/gj- belongs toiP{iF). In particular, 

best n-term Hermite approximations converge in L^(C/, C, 7 ) with rate n~^ where s = | ~ 5 = |- 

As already observed, we always have q{p) > p, which shows that there is in this case an im¬ 
provement in the summability properties of (||u,y||\/)j/gjr over those of (||V’j||L°°)j>i- For example, 
summability of (||V'illl/°°)i>i implies summability of (||ui/||y),yg:F j and therefore convergence 
of best n-term Hermite approximations with rate 


6.2 Wavelets 

For a general wavelet bases on a domain D we adopt the notation {'fix); used for example in 

[ 6 ], where A concatenates the scale and spatial indices, with the convention that the scale level I of 
■fix is denoted by |A1, i.e., |A| := 1. Thus, there are 0{2‘^^) wavelets at level I and the wavelets at 
each given scale have finite overlap. This means that for all x € D, 


ff{X : |A| = / and fix{x) 7 ^ 0} < M, 

for some fixed M > 0 independent of 1. We consider wavelets normalized such that 

,^=ci = C 2 -“^ |A| = I, 


(127) 


(128) 


for some fixed C > 0 and a > 0. Using the finite overlapping property, we find that for any 
0 < K < a the sequence 

Pa:=2"I^I, (129) 


satisfies 


sup ^ Px\Mx)\ <CMJ2 < 00 , (130) 

xGD 


A l>0 

which is (25) in Theorem 1.2. 

Note that if we order our wavelet basis from coarse to fine scale, we find for the resulting system 
{'fij)j>i and sequence {pj)j>i the algebraic behaviour 


bj\\L'=° ~ j 


— OLjd 


(131) 
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and 


( 132 ) 


We thus obtain from Theorem 1.2 the following immediate corollary. 

Corollary 6.2 Let {'4’j)j>i be a wavelet basis with the normalization (128). (HV'j||L°°)j>i belongs 
to then (||rti,|| belongs to £P(iF) for allp such that > q. In particular, the best n-term 

Hermite approximations converge in L‘^{U, V,'y) with rate for all s < |. 

Proof: Since {\\'i(j\\L°°)j>i belongs to we know that qa > d. Since for the given p, we have 

q* := > q, we can take k and {p\) as defined in (129) for this k. Then, the conditions 

of Theorem 1.2 are satisfied for this sequence of weights and for q*. Hence, (||n,y||y)i/gj- belongs to 
LP{IF). This also gives that the approximation rate is with s = 1/p — 1/2 = I/q*■ By adjusting 
p, we can take q* < q as close to q as we wish, thereby establishing the rate of convergence n~^ for 
all s < l/g. □ 

Note that if we use sufficiently smooth wavelets, the decay property (128) is equivalent to the 
property that the correlation function belongs to the Besov space B^{L°°{D)), which coincides 
with the Holder space C" when a is non-integer, and therefore b is almost surely in the Holder 
space for /3 < a/2. Thus, we also infer from Theorem 1.2 that if Ch belongs to the Besov space 
for some a > 0, best n-term Hermite approximations converge in L‘^{U, V,'y) with rate 
n~^ for all s < a/d. 

6.3 Arbitrary supports 

We finally consider functions ipj with arbitrary supports, including the case of globally supported 
functions such as the Fourier basis. Let us assume that (H'j/i l|L°°)i>i € for some 0 < g < 1. 

We then find that the choice 

Pj := llV’jllioo , (133) 

obviously satisfies (25). We also find that {pj^)j>i belongs to £'’(N) for r := Therefore, 

applying Theorem 1.2, we find that (||njy||y)i/gj- £ when p satisfies r = or equivalently 

I = 1 Therefore, we obtain the following immediate corollary. 

Corollary 6.3 Let (V’i)j>i be a family of functions with arbitrary support, and let 0 < p < 2. If 
(||V'illL°°)i>i belongs to with | = 1 + ^, then (||ni/||y)i/gjr belongs to In particular, best 

n-term Hermite approximations converge in L^([/, H, 7 ) with rate for s = | — 1. 

Note that if (j^||V’illl/°°)i>i £ ^^(N) for some /3 > ^, an application of Holder’s inequality shows 
that (llV’j|| l°°)j>i belongs to £'^(N) with | = ^ + ^- Therefore, the above corollary represents an 
improvement over the condition (iUV’jl|L°°)i>i ^ ^^(N) from [19]. 
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7 Non-optimality of the Karhunen-Loeve representation 

The previous examples illustrate the role of the support properties of the functions {'4’j)j>i when 
analyzing the convergence rate of best n-term Hermite approximations for the map y i—)■ u{y). In 
particular, they reveal that faster convergence rates can be obtained in the case of locally supported 
functions. 

We now discuss a concrete example which illustrates this phenomenon for the approximation 
of a PDE with given lognormal coefficients, when we use two different representations of these 
coefficients. Here we take 

D=]0,1[, (134) 

and therefore consider the equation 

— {au'Y = /, u(0) = tt(l) = 0. (135) 


We take a = exp(6) where 6 is a Brownian bridge, that is, a Gaussian process with covariance given 

by 

Cb{x,x) = min{x,x^} — xx . (136) 

There exists two simple explicit representations for this process. 


1. The Karhunen-Loeve representation is determined by the eigenfunctions and eigenvalues of 
the covariance operator, which have the form 


(Pj{x) := \/2 sin(7rjx) and Xj := ^ , j > 1, 

so that after normalization, we obtain the representation 

a/2 

Hy) = ^ yjY>j: i’jix) ■= — sin(7rjx), 
i>i 


(137) 


(138) 


where the yj are i.i.d. standard Gaussian random variables. 

2. The Levy-Giesielki representation uses the Schauder basis functions, which are the primitives 
of the Haar functions, namely 

^i,fc(x) = 2-'/2^(2 'x -/c), A: = 0,...,2'-l, Z>0, (139) 

where V’(^) := rnaxjO, 1/2— \x— 1/2|}. Then setting ijjj = when j = 2^ +k, we have again 

i>i 

where the yj are i.i.d. standard Gaussian random variables. 
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If we use the Karhunen-Loeve representation, the analysis carried in this paper does not allow 
us to establish any summablity of for p < 2. Indeed, due to the global nature of 

the functions we need to rely on Corollary 6.3 which requires that i\\'4’j\\L°°)j>i £ ^^(N) for 
some q < I, which does not hold since this sequence is not summable. Note however that 
summability is ensured since E(||u(y)||y) < oo. 

In contrast, when using the Levy-Ciesielski representation, which is essentially of wavelet type 
with a = ^ in (128), we can rely on Corollary 6.2 which shows that belongs to 

for all 1 < p < 2. In particular, best n-term Hermite approximations converge in L^(17, fo, 7) with 
rate n~^ for all 0 < s < 1. 

This example reveals that for a given lognormal process, the Karhunen-Loeve representation of 
the Gaussian process might not be optimal in terms of the resulting convergence rates of the best 
n-term Hermite approximation. One heuristic explanation of this fact is that the Karhunen-Loeve 
representation is optimal in a very specific sense: it minimizes the mean-square L^(H)-error when 
truncating b by the J first terms in its expansion. However, in the present setting of the elliptic 
diffusion equation, the relevant norm for approximating the functions a and b is not the norm, 
but rather the L°° norm for which the Karhunen-Loeve representation has no particular optimality 
property. 
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